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O ' 

^— ^ . Continuing our development of a mathematical theory of stochastic microlensing, we study the 

^sj ' random shear and expected number of random lensed images of different types. In particular, we 

characterize the first three leading terms in the asymptotic expression of the joint probability density 

1^ . function (p.d.f.) of the random shear tensor due to point masses in the limit of an infinite number of 

stars. Up to this order, the p.d.f. depends on the magnitude of the shear tensor, the optical depth, 
and the mean number of stars through a combination of radial position and the stars' masses. As 
a consequence, the p.d.f.s of the shear components are seen to converge, in the limit of an infinite 
number of stars, to shifted Cauchy distributions, which shows that the shear components have heavy 
tails in that limit. The asymptotic p.d.f. of the shear magnitude in the limit of an infinite number of 

r^ , stars is also presented. All the results on the random microlensing shear are given for a general point 

f~^ I' in the lens plane. Extending to general random distributions (not necessarily uniform) of the lenses, 

we employ the Kac-Rice formula and Morse theory to deduce general formulas for the expected 
total number of images and the expected number of saddle images. We further generalize these 
results by considering random sources defined on a countable compact covering of the light source 
plane. This is done to introduce the notion of global expected number of positive parity images due 
to a general lensing map. Applying the result to microlensing, we calculate the asymptotic global 
expected number of minimum images in the limit of an infinite number of stars, where the stars 
Cn ' are uniformly distributed. This global expectation is bounded, while the global expected number 

^ ' of images and the global expected number of saddle images diverge as the order of the number of 

\l , stars. 

00 
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I. INTRODUCTION 

We started an analytical study of stochastic microlensing in Paper I (Fetters, Rider, and Teguia [l|), where we 
^ ' derived stochastic properties of the random time delay function, lensing map, and bending angle. We continue this 
work here by investigating the expected number of critical points (lensed images) of the time delay function, or 
equivalently, the expected number of solutions of the random lens equation. 

There is a fair amount of literature dealing with the non-random number of images produced by a microlensing 
situation. One of the first works on this subject showed that a fixed source lensed by a single star can produce two 
distinct images (Einstein [4]). This result was generalized to any finite number of stars for both single and multiple 
plane lensing, with results on sharp lower bounds and counting formulas of the total number of images, the number 
of minimum images, number of saddles images and number of maximum images (Fetters [3|, |j, |5[; also see Chapters 
11 and 12 of the monograph by Fetters, Levine, and Wambsganss [6j). Sharp upper bounds on the total number of 
images have also been obtained (see S. H. Rhie [7] and Khavinson and Neumann [8|). The case of stars with random 
positions was later examined using semi-analytical methods, yielding asymptotic formulas in the limit of an infinite 
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number of stars for the expected number of minima (e.g., Wambsganss, Watt, and Schneider y\; Granot, Schrecter, 
and Wambsganss [i3|)- However, these results are hmited to first order approximations and are not appUcable to 
general random distributions of the star positions or source positions. 

In the current paper, we begin by deriving the first three terms in the asymptotic expansion of the microlensing 
shear tensor's p.d.f. in the limit of an infinite number of stars. These results show that in such a limit, the marginal 
p.d.f.s of the components of the random shear tensor converge in the sense of distributions to shifted Cauchy densities. 
We also show that the third term in this asymptotic p.d.f. of the shear tensor depends on the square of the magnitude 
of the shear due only to stars, the optical depth, and the mean number of stars. The asymptotic p.d.f. of the 
magnitude of shear is determined as well. 

Next, we derive for a lensing map (not necessarily due to microlensing), a Kac-Rice type formula for the expected 
number of positive parity lenscd images (i.e., minimum and m.axim,um, images) lying in a compact region of the 
lens plane. This formula can be applied in a large number of random lensing scenarios, and is obtained using the 
Kac-Rice technique for computing the expected number of zeros of random mappings (e.g., Adler and Taylor [ll|; 
Azais and Wschebor [l2|; Forrester and Honner [ij]; Sodin and Tsirelson [3, [la [3; Shepp and Vanderbei [l3])- We 
then combine this formula and results from Morse theory to obtain additional general formulas for the expected total 
number of images and the expected number of saddle images. 

The paper then concludes by applying the aforementioned Kac-Rice formula and the asymptotic joint p.d.f. of the 
microlensing shear tensor to microlensing due to point masses with continuous matter and shear. Here, we focus our 
attention on the "global expectation" of the number of minimum images (which will be defined later) . We consider the 
minimum images since numerical simulations have shown that they are significant contributors to high magnification 
events in microlensing light curves (e.g., Wambsganss, Watt, and Schneider Q). 

The outline of the paper is as follows: The basic lensing framework is given in Section |TT1 We obtain the asymptotic 
joint p.d.f for the microlensing shear components in Section IHII In Section llVi we derive general formulas for the 
expected total number of positive parity images, the expected total number of images, and the expected total number 
of saddle images. In the same section, we define the global expectation and obtain an integral formula for the global 
expected number of positive parity images. Section |V] presents applications of these general formulas, along with the 
results on the random microlensing shear, to determine a global expected number of minimum images in stochastic 
microlensing. 

II. BASICS 

A general gravitational lens system is governed by its time delay function 

and lensing map 

ri(x) — X — Vf/'(a;), 

where y is the light source's position and ij} is the gravitational lens potential (Ref. [6|, Chapter 6). A lensed im,age of 
a light source at y is a critical point x oi Ty, that is, a solution x in R^ of the equation VTy{x) = 0. It immediately 
follows that lensed images are solutions x in R^ of the equation 

r]{x) = y, 

which is called the lens equation. When the potential ip or the source position y becomes random, the lens equation 
has a random number of solutions (lensed images) . 

A lensed image x is nondegenerate (resp., degenerate) if a; is a nondegenerate (resp., degenerate) critical point of 
Ty. Using the Morse lemma (see Milnor |l31), we can classify a nondegenerate light ray as either a local minimum (if 
X has index 0), a local saddle (if x has index 1), or a local maximum (if x has index 2) (see Ref. |6j, p. 178 for more 
information). The set of positive parity lensed images consists of all the minimum and maximum lensed images. 

In stochastic microlensing, we consider a gravitational field due to a collection of g stars at positions fi, . . . , ^g that 
are independent and identically distributed (i.i.d.) over a region of the lens plane. The resulting gravitational lens 



potential f/jg, time delay function Ty_g at y, and lensing map 77^ become random and are given by: 



V'g(x) 



-{xl-xl) + ^ rrij log \x - ^^ 



j=i 



1 • ® 



5 



%(2;) = ((1 - Kc+ l)xi, (1 - Kg + 7)xi) - y^ i 



J = l 






where a; = (a;i,a;2) is in L = M^ — {^1, . . . ,^g} (lens plane). Unless stated to the contrary, the light source positions 
y, the continuous matter density Kc, external shear 7, masses rrij, and number of stars g are assumed fixed. The star 
positions ^j = {Uj, Vj) are random vectors. We shall consider the case of random light source positions in Section HVl 
and Section IVl 

In Section [ml and Section |Vl we shall work under the following assumptions about microlensing: 

• Cartesian coordinates {u, v) are assumed in the plane M^ with as the origin. 



• All masses are assumed equal: mi — 1112 



m„ 



m ^ 0. 



• The star positions ^i — {Ui,Vi) are independent and uniformly distributed in the disk B{0,R) of radius R 
centered at 0, that is, 



{U^,V,)^Vmi{(u,v) gK^ : \u\^ + \v\^ < R^} , « = 1, 



,9- 



• The quantities R and g are related by the following physical formula for the optical depth k* in point masses 

(stars): 



^^ — 



m g 
i?2 • 



(1) 



The quantity k* is fixed. For x E R^, with < |x| < R, we shall call the combination K*|a;p/m the mean 
number of stars within the disc of radius |x| about the origin. 

Remarks: The assumptions stated above provide a good approximation for several astrophysical scenarios. Note 
that another lensing model used in the Hterature is the elliptical lens model (see Ref. Q, PP- 101-108 for more on 
lens models). 



III. RANDOM SHEAR DUE TO POINT MASSES 



The shear at x is a gravitational tug of stars on the light ray through x. The two components of the shear tensor 
at X = (xi,a;2) G L due to point masses (stars) is denoted by T*{x) = (r| (x),r2 „{x)), with components as follows 
(see Ref. (I, p. 96): 



^1,0(2^) — y. Trrr \2~~7T7 \2i2~ ^"^^ ^2,g{^) 



3 



2m{Uj — xi){Vj — X2) 

^^[{u,-xir + {v,-x2w 



We study the asymptotic behavior of the random shear tensor {Tl „,T2 „) due to the point masses by obtaining 
several leading terms of the asymptotic joint p.d.f. of its components in the limit g —^ 00. To do so, we assume (for 
the remainder of the paper) that the random shear tensor has an absolutely continuous distribution function. We 
then obtain the following: 

Theorem 1. Let x — (a::i,a;2) G M^ with \x\ < R be fixed. The joint p.d.f. /p* (x).r' (x) of the pair {Tl „(a::),r2 n{x)) 
is given in the limit g — + 00 by: 



/r* (x),r; (x){z,w) = „ , „ , t^ox^/o 



i?i(r) ^ H2{T; \x\) 



+ Oig-'), 



(2) 



where {z,w) denotes the possible values of (F^ q{^)t^2 qi^))' T — \J z^ + lip- denotes the possible values of the shear 
tensor's magnitude, and 



Hi{z,w) = Hi{r) = k: 



2 9r"-6Ac; 

4(k2 + r2)2 



jy2(2,«;;la:|) = HiiT; \x\ 



K,\xf ^2(6^2 _ gp2-, ^2(g^4 _ 24^2p2 ^ 3p4>, ^5^4 (^g^4 _ 4Q^2p2 _^ ^5p4 



m 2(«:2+r2)2 



4(«2 + r2 



+ 



32(^2 +r2)4 



(3) 



(4) 



The p.d.f. at this order is a function of the magnitude of the shear due only to stars, optical depth k^, and the mean 
number of stars within the disc of radius \x\ about the origin. 
The marginals of f^* {x),r' i^) ^'"^ 



/r,,g(2;)(ai) 



7r{Kl+af) 



1 



Xi(a2) , K2ial, \x\) 



+ 0(5-3), i = l,2, 



where ai — z, a2 — w and 

Kiih) = kI 
K2ih; \x\) = 



3h-Kl 
2{Kl + h)^- 
«;*|a;pK*(K*-3/i) KlKt-6Klh + h^ Sut 4 - IOkIH + bh"^ 



m {kI + hy 



(k2 + hf 



{nl + hY 



Proof: We set up coordinates in the same way as in the proof of Proposition 1 in [l|. For the reader's convenience, 
we repeat the details. Let x = (a;i,a;2) S M^ with |j;| < R. Translate the rectangular coordinates {u,v) so that 
its origin moves to position x and denote the resulting new coordinates by {u',v') and its origin by . The old 
origin in {u' ,v') has coordinates {u',v') = —x. Define loq to be the unique principal angle with coscjq = a;i/|a;|, 
sinwo = a;2/|a;|, and let w = wq + tt. Rotating the rectangular coordinates (u',v') counterclockwise by angle lu we 
obtain new coordinates {u",v") with origin 0" = 0'. The old origin now lies at position {u",v") — {\x\,0) = x" on 
the positive M"-axis. Let {9,r) denote polar coordinates in the frame {u",v"). 

Let X — (xi,X2) € M^ with |a;| < R and set Ux ^ u — xi, v^ = v — X2, and r^ = m^ + w^. The joint characteristic 
function of (FJ_g(a;), F2,g(x)) is given by: 



{'Pri^{x).rijx)iti,t2) 



1/9 



7ri?2 



7ri?2 



exp 



B(0,B.} 



■ . , m{ul-vl) \ . .. (2rriu^ 



dudv 



exp 



B{x,R-\x\) 
1 



.m 



i^{ti,t2) ■ iux-v^,2u^v^) 



dudv 



7ri?2 



exp 



i?2 



B{0,R)\B(x,R-\x\) 






rJo 



dr 



7ri?2 



exp 



imt cos {2(9 + to')) 



dudv. 



dOdr, 



(5) 



with t = ^Jtl+tl, 2cc;' = 2lo - 
Now, with V = mt, we get: 



— arctan |^ , Jo is the zero Bessel function. 



rJ( 



(^) dr 



v 

i?2 



1^2(.--,-,i 



U 



2' 2' 'A{R-\x\ 



v_ {R-\x\f r 

i?2 ^ i?2 



4(i?-|.T|)4 



0{R' 



(6) 



where 1-F2 is the hypergeometric function. 



We use Taylor series expansion and term-by-term integration to write the second integral in equation ([5]) as: 



1 



cxp 



iiy cos {2{9 + Lu')) 



where, for integers n > 1, 

r-fl+bl rfir) cos" (2(0 + w')) 



dudv = 



7ri?2 



7:{2\x\R - l^n + ^l!^/„ 



n=l 



In = 



R~\x\ J-f(r) r 

R+\x\ 

R-x 



2n 



dedr 



1 /■^('■^ f ^ \ 

— T / «« + V c„j- cos [2(n - 2j)(0 + u')] dedr 



R+\x\ , 

:5;^|2a„/(r)+^ 



Cn.j COS [2(n — 2j)w'] . 



c„.o COS [2naj'] 



.,, -- L^v- - ^~ J. sin [2(n - 2j)/(^)] I dr + -"■» — ^— ^ 3^^ 



i?-k 



n-2j 



where 



R-\x\ ' 

and the coefficients a„, 5„ and c^^jS are obtained by finding the Fourier series of cos" 0. We will prove in the appendix 
that 3„ = for every integer n > 2. We then obtain: 

h - 0, 



R-\x\ 



2 (i?2 „ |2:|2) 



2^2 



^^ ^ 3cos2co^ |-^+l"l cos(/(r)) sin(/(r)) ^^ ^ 3cos2co^ 7r\x\^R^ ^ 0(R-') 



R-\x\ 



2 2(i?2 _ |2;|2) 



2M 



|/„l < 



R+\x\ 



R-\x\ 



'^ + 'dr < ^'!'(^\+\) = Ofi?-"+^) forn> 0. 



„2n-l 



(i?-|a;|)2»-i 



Therefore, 



7ri?2 



exp 



^008(2(6* + ^;')) 



dudv 



2\x\ \x\^ \xW' \xYv' Wv 



2, ,2 1^13, ,2 



R i?2 2i?5 4i?6 ^7 



o(i?- 



(7) 



We now substitute equations ([6]) and ([7]) into equation (O to obtain 



\fTl^^{x),Vl^^{x){txM)) 



^/g V V 



|.t|2z/2 



1 



i?2 4i?4 2i?6 

K,t K2f2 \x?-Klt^ 



g Ag^ 2mg'^ 



OiR-"^) 
+ 0{g-') 



An asymptotic expansion of the characteristic function is then derived apply the logarithmic function of both sides 
of the above equation: 



log(^^rt,,(a:),r5,,(x)(^i,i2)] = 5 log 



1 + -V + -4 — V + 0(g ) 



g Ag^ 2mg^ 



thus. 



''Prijx)sijx){h,t2) = exp 



4,9 



|a;|2K*j/2 jy3\ 2 



= 1- ^ 
45 



2™ -r2)7^''^'\ 

2m 12 32/ g^ 



Hence 



VTl,{x).Tl^(x){tl:t2) - e" 



4g 



\xYk.-j,vi Vi 



2m 



12 



32/ g 



0{g- 



(8) 



By assumption, the joint p.d.f. /p* (a;),r' (x) of ^i 0(2^): ^2 q{x)) exists, i.e., the Inverse Fourier Transform (I.F.T.) 
of the L.H.S. of equation ([5]) exists. Consequently, so does the I.F.T. of the R.H.S. of equation ([5]). The later can be 
computed as the sum of the I.F.T. of functions of the form p{g) x (?(ti, 12) where p{g) is a function of g and q{ti, ^2) 
is independent of g and integrable. 

Let 



'Pg{h,t2) = er 



1 



43 



\x^K^v\ v\ 



1m 



12 327 52 



1 



We have 

/ri„W,r5„w(^'H = (IFT(^<,))(z,u;) + 0(5-3). 
Here IFT((^g) is the I.F.T. of ipg and can be obtained by computing integrals of the form 

1 



(2^)2 



27r 
1 



(z2+u;2)~n/2-l / ^n+l^-K. (-^+-^)-^'^ -Jq (r)dr 



We then complete the proof of the first part of the theorem by using equation ([9 
Now consider the following functions defined on M: 

h\{z) = fri^(x)x^^{x){z,w) for fixed weR. 
h2{w) = fri^(x)xi^ix){z,w) for fixed z e M. 



z2 + w^ 



for n = 0,2, 3, 4. 



kliz) E 


K, 1 

" 2^(k2+z2+w2)3/2 


'_! Hi{z,w) H2{z,w; 
[ 5 ' 5^ 


\A) 


k2iw) E 


K, 1 

27r(K2+22_^y;2)3/2 


r 1 Hi(z,uO 1 1/2(2, tf^; 
5 9^ 


N) 



(9) 



for fixed w G 
for fixed z G I 



We know that hi, ft.2 G i^(R) since the random shear's p.d.f. is integrable. We can then prove that fci, ^2 G i^(K) 
by direct computation. Hence, the second part of the theorem follows by integration with respect to w and z. D 



Remark: The leading function in the R.H.S. of equation ([2]) becomes a p.d.f. asymptotically in g. In particular, set 



^gAx\{z,w) 



Then for g sufficiently large. 



In addition, 



Hi{z,w) H2{z,w; \x\ 



2tt{kI + z^ + w^)y^ 

^g,\x\{z,w) > 0. 

^g, \x\{z, w)dzdw 



for all 5 > 1 and a; G M^. An analogous result was found in [l| for the random lensing map (see the discussion after 
Corollary 10). 



Corollary 2. For any x G M^, the pair {r\ Ax) ,T2 „{^)) converges in the sense of distributions to the pair 
(rj oo,r2 oo) with joint density given by: 



/] 



rr ,.,r* 



(z,w) 



2tt [kI +z2 + w2]3/2 
as 5 —!■ oo. Note that each marginal of fy ,y* is a {stretched) Cauchy distribution: 



friJ^) 



Tt{kI +z2) 



and fr- ^ (w) 



7t{kI + Up- 



Corollary [5] is consistent with the findings by Nityananda and Ostriker [10| and Schneider, Ehlers, and Falco [201 
(see p. 329). Note that the corollary shows that the random shear components due to point masses have heavy tails 
as 5 — > cx), an important point we have not seen emphasized in the literature. 



Corollary 3. Let a; G K^ with \x\ < R be fixed and consider the magnitude of the shear at x due only to stars, namely, 
^9^^) = V (^i.s^^)) + (ri,g(^)) ■ ^^^ P-'^-f- fr;{x) ofT*g{x) is given in the limit g -^ oo by: 

, Hi{T) , H^iT; \x\ 



/r-(.)(r) 



/^*1 



(K2+r2)3/2 

Here T denotes the possible values ofT*{x). 



+ 0(9-^). 



Proof: We use the change of variable 



and obtain 



Tljx) = r;(x)cose(a:) 

n,gi^) = r;(x)sine(x), 



fr'(x),e{x){^,S) = r/p. (^),r. (:r)(rcos6',rsir 



2^(K2+r2)3/2 



gi(r) H2ir;\x\) 

9 

9 9 



+ 0{g-'). 



We can then integrate the equation above with respect to 6 from to 27r and obtain the result stated in the Corollary. 

D 



IV. EXPECTED NUMBER OF IMAGES: GENERAL CASE 

In the current section, we are not restricting ourselves to the point-mass model as was the case in the previous 
section. We shall determine results for the expected number of images for a general lensing map. 

A. Expected Number of Positive Parity Images: General Case and the Kac-Rice Formula 

We work under the following broad, natural assumptions and notation: 

• The gravitational field potential ip is an isolated potential with a finite number of singularities, all of which are 
infinite (see [g], p. 424 for more on isolated potentials). 

• The potential tp and lensing map 77 are smooth, except on a set of measure zero. 

• Notation: 

' 1 if X e A 



For a set A, let 1a{x) = ^ ^ if x (^ A. 

— Set G{x) = det(Jac?7)(x), A = (0, 00), and tj^ = {i^ G M^ : G(i^) G A}, where i] a general lensing map. 

— Let iV+ = N„iin + Nmax, which is the total number of maximum and minimum images. 

— Let N^(D, y) — total number of positive parity images (minima and maxima) lying in a subset D of the 
lens plane for fixed source position y. 

The determinant of the Jacobian of the general lensing map 77 is given by (see [6| , p. 182): 

det(Jac?7)(a;) = (1 - n{x)f - T^-{x) = (1 - K{x)f - [tI{x) + TI{x)\ . 

Here Fi and F2 are the shear components of the general lensing map 77 and are defined by 

Fi(a;) ^ and F2(a;) = -012(3;), 

and K,(x) is the convergence defined by 

i)\\{x) -|-i/'22(a;) 



k{x) 



where i/jy — d'^'^jdxidxj. 
The main result is: 



Theorem 4. (Expected Number of Positive Parity Images) Let D cM? be a closed disk. Assume the following: 

1. P[{xeD : ri{x) = y and det(Jac77)(2:) = 0}] = 0. 

2. K is bounded on L {lens plane). 

Then for a test function p :M.^ —^ W with compact support, we have 



p{y)E[N+{D,y)]dy ^ / p{y) E dei{Jacri){x) lg^{x) ri{x) ^ y f^^,){y) dx dy. (10) 

and, for almost all y m R^, the following Kac-Rice formula holds for lensing: 

E[N+{D,y)]= [ EU{l-n{x)f-V\x)) lg,(a;) ii{x) = y] fr,(.){y) dx. (11) 



Before we prove this theorem, it is important to note that it can be applied to a large range of lens models. 
Moreover, the singularities form a set of measure zero and so contribute a value of zero to these integrals (see Rudin 

im, p. 20). 

The first assumption of Theorem 3] requires that the random images are almost surely not critical points of the 
random lensing map. The second assumption holds for optical depth constant in L. 



Proof: The proof follows, with some modifications, the proof of the Kac-Rice formula given by Azais and Wschebor 
in [12] (Theorem 2.1). A more general version of the Kac-Rice formula with a different proof is given by Adler and 
Taylor in [ll|. The main tool used in the proof below is Federer's co-area formula. 

We start by fixing the source position at y. By assumption 1 and results in [3, |j, 0, Q, the random function 
N+{D, y) is uniformly bounded, except on sets in M^ of measure zero. Let p : E^ ^ K be a test function with compact 
support. Since rj is smooth almost everywhere in M^ and p{y) is integrable, we can apply Federer's co-area formula 
(see Federer [ll|, p. 303) and obtain 



p{y)N+{D,y)dy^ / \det{3acri){x)\p{r]{x))lg^{x)dx. 
Jd 

Taking expectations on both sides of the above equation yields 

piy)E[N+{D,y)]dy= f E[dct{3a.cr^){x) p{ri{x))lg^{x)]dx. (12) 

The L.H.S. (and hence the R.H.S.) of equation (fT2)) is bounded and the absolute value in the R.H.S. is dropped since 
we are only considering positive parity images. Using the basic result 



E[U] = f i?[t/|V = v]/v(v)dv, 



(13) 



where [/ is a real- valued random variable, V is a M"-valued random variable, and v G M", we obtain: 



£;[det(Jac?7)(2:) p{r]{x))lg^{x)] = 



E 



det(Jac77)(a;) p{Tj{x))lg^{x) r](x) ^ y f^(x){y)dy 



Note that by equation p^ . the expectation i? [det(Jac77)(a;) p{ri{x))lg^{x)] is bounded on D, except on subsets of 
measure zero. The above equation then implies that 



D JR2 



E 



det(Jac77)(a;) p{i]{x))lg^{x) ■q{x) ^ y f,^(^){y) 



dydx 



is finite. Therefore: 

p{y)E[N+{D,y)]dy = 



2 Jd 



D 



det(Jac77)(a;) p{T]{x))lg^{x) t]{x) = y f^(^){y)dydx. 
det(Jacr7)(a::) p{'q{x))lg^{x) ■q{x) = y fn(x){y)dxdy 
det(Jac77)(x) p{y)lg^{x) r]{x) = y fr,(x){y)dxdy 



p{y) / E det(Jac77)(a;) l£;^(x) r]{x) = y .fr,(x){y) dx dy, 
Jd l J 

where the second equality follows from Fubini's theorem, the third follows from the condition ri{x) — y, and the last 
one follows from the fact that p{y) is deterministic and independent of x. This proves the result in equation (|10p . 
Since the test function p was arbitrary, equation pTjl is valid for almost all y e M^. D 



Remark: For the critical case, if k(x) — 1 for almost all (a. a.) x E L, then G(x) ^ A for a. a. x G L, so lg^(a;) = 
for a. a. x € L. Hence, E{N+{D,y)) = 0. This is consistent with results in Ref. Q. 



B. Expected Number of Images and Saddle Images: General Case 

We can modify Theorem 2] to obtain a formula for the expected total number of images and the expected number of 
saddle images lying in D. However, it is simpler to use results from Morse theory contained in the following theorem 
(see [1,13 and (H], Theorem 11.4, p. 424): 
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Theorem 5. 6^ Let T and 77 he, respectively, a single-plane time delay family and lensing map induced by an isolated 
gravitational lens potential 7p. Suppose that ip has g singularities, all of which are infinite singularities. Then the total 
number of lensed images of a light source at a non-caustic point y is finite and the following holds: 

1-N^,n>l, Nsad>9 + N^ax, N^ad > N^ir, + g - I, N+ = N^ad - 9 + I- 

2. Ntot = 2iV+ + g - 1 = 2Nsad - ff + 1, Nt,t>9 + 1- 



For a generalization of Theorem[5]to the case when the potential is not isolated, see Ref. [6|], p. 422. We can apply 
Theorems [S] and m to obtain formulas for the expected total number of images and expected number of saddle images: 



Corollary 6. Let D be a closed disk whose interior contains the g singularities ofip. Then, under the assumptions 
of Theorem pi the following holds for almost all y in R^ ; 

E[Ntot{D,y)] = 2 / E\{il-K{x)f-T\x)) Igjx) v{x) = y] f^^M dx + g ^ 1 



D 



E[NsadiD,y)] 
with self-evident notation. 



E 



D 



{{I- K{x)f ~T'^{x)) Ig^ix) ■q{x)^y f^(^){y) dx + g 



Proof: Let ^i , • • • , ^g be the positions of the singularities. Choose eg small enough so that no image lies inside B{^j , eg) 
for j = 1, . . . ,g and B{£^j, cq) O B {£^i , eq) =0 for all i y^ j. This is possible since the total number of images is bounded 

by 

g + l<Ntot< 5.9 - 5 

(see [3, Ijl for the lower bound and 0, Q for the upper bound) . 

Set D' = D\[J^--^^B{^j,€o). The time delay function Ty associated with V' generically satisfies Morse boundary 
condition A relative to D' (see , p. 403) . We can then apply Theorem [5] to D' and obtain 

iVL = 27V; + .9 - 1 = 2iV;,^ -5 + 1, 



(14) 



where quantities with "prime" are related to D' . 

Equation (J14p is an equality of random variables, and holds for almost all sources position y e R^. Therefore, we 
can take the expectations and obtain: 

E[Ni,] = 2E[N'^] +g-l = 2E[NU -g + 1. 

However, by construction, 

Ntot{D,y) ^ Ni„ N+{D,y)=N'^, and iVw(Ay)=A^w- 

The results then follow from Theorem [H D 



Remark: In Corollary [6l the expected total number of images and number of saddles images both depend on the 
topology of the region DHL since g is its first Betti number. 
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C. Global Expectation: General Case 

The expectations we shall study are actually generalizations of the expression E[N+{D,y)] in equation pT|) since 
the source position is not known in the physical systems of interest to us (e.g., microlensing) . Intuitively speaking, 
the quantity E\N j^{D , y)] incorporates averaging with respect to the random components of the lensing map r] for a 
fixed source positon y. We shall extend this to averaging over random source positions in a fixed compact set with 
positive measure in the light source plane and then generalize further to "averaging" over a family of such compact 
sets with positive measure. When all these averages are carried out, we refer to the result as a "global expectation." 
We now make these ideas precise. 

The first extension of E[N^{D,y)] is to consider a light source position in M^ that is not fixed at y, but whose 
position is a random vector Y with density fy ■ We then define the expected number of positive parity images in D 
as follows: 

E[N+{D,Y)]= f E[N+{D,y)] fY{y)dy. (15) 

When the p.d.f. fy of Y has compact support &o, we write the p.d.f. as /y^eo to emphasize the support and denote 
the expectation (fTS]) by: 



E[N+{D,Y:eo)]^ [ E[N+{D,y)] fY,eMdy. 

JM2 



We have: 



Corollary 7. Let Y be the random source position having a general p.d.f. /y,eo vjith compact support ©q. Then, 
under the assumptions of TheoremlA we obtain: 

E[N+{D,Y; 6o)]= f E [{{1 - ^ixjf -T^x)) Igjx) fy^eoivi^))] dx. (16) 

Jd 



Proof: The result follows from equation P^ with p — fY,eo- D 



Remark: It is important to add that the general formula in p6|) does not require knowing the distribution of the 
pair (ri,r2) conditioned on r/ = (771,7^2), as is the case for equation (jlip . But equation (|16p requires knowledge of 
the joint distribution of the quadruple (Fi, T2,r]i, 772)- However, as will be seen in Section |Vl knowledge of the joint 
distribution of the pair (Fi, F2) will be enough for a certain microlensing case of interest to us. 

Example. Suppose that Y is uniformly distributed with p.d.f. fY,6oy where the compact support 60 is a subset of 
the light source plane with nontrivial area: 



Then: 



area(6o) = |6o| > 0. 



E[N+{D,Y;e,)] ^T^ E [((1 - K{x)f - T^{x)) lg,(a;) leM^))] dx. 
|t'o| Jd 



We are now ready to state the generalization of expectation from the set ©o with positive measure to a family of 
sets. Let {©} be a countable family of compact subsets with positive measure covering the light source plane M^. 
Suppose that each member © has the same area as a fixed member ©q. For each 6, consider random source positions 
Y that are uniformly distributed with p.d.f. fy^e compactly supported in ©. Note that by construction we have: 

r I X i6(y) i6(y) 
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Now, the quantity E[N+{D, Y; &)] captures an average with respect to the randomness arising from the lensing map 
rj as well as an average with respect to the random source position Y over 6. We now "average" E[N+{D, Y; &)] over 
the family {&}. Explicitly, we define the global expected number of positive parity images lying in D for a uniform 
source by: 

E[N+iD,Y;e)][e} = T^ / S [det (JacT?) (x) la^(x)] dx, (17) 

Ic'ol Jd 

where \6\ = |So|- The global expected number of minimum images E[N^{D,Y; 6)]{e} is well-defined (i.e., bounded) 
since 

< det ( Jac r]) (x) Ig^ (x) < 1 

and D is compact. This definition captures the physical setup we need for the microlensing application in Section IVl 
We shall compute (|17p asymptotically to three orders in the case of microlensing. 

V. GLOBAL EXPECTED NUMBER OF IMAGES: MICROLENSING CASE 

We now consider the case of microlensing, which consists of randomly distributed point masses with continuous 
matter and shear, as described in Section |TT1 We shall apply the results given in Sections [Til and [TVl to study the global 
expected number of minimum images in microlensing. 

Consider a closed disk D. The random lensing map rig satisfies the assumptions made at the start of the previous 
section and in the statement of Theorem [H In fact, 

P[{x £ D : rjg{x) = y and det(Jac77g)(a;) = 0}] = 0, 

(i.e., lensed images in D are almost surely not images of caustic points) and 

P[dD C\{xeD : f]{x) == y} = 0] = 1 

(i.e., there is almost surely no image on the boundary of D) since there arc finitely many lensed images for almost all 
source positions y (see p], p. 431 for more details). Note that for microlensing, we have 

ri{x) =^ "/ + Tlg{x), r2(a::) = r2_^(a;), and k{x) ^ Kc ior x € L, 

where L — M."^ — {^i, . . . , ^^j (lens plane), and that D may contain singularities £,i. 

Assume that Kc < 1 (subcritical) in the current section. This implies that Nmax = (see [y], p. 433). Consequently, 
for the microlensing map rjg with Kc < 1, we have the following for any non-caustic source position y: 

N+{D,y) = N,n,n{D,y) = Ng,^,^{D,y). 

Furthermore, the set tj^ = {i^ € R^ : det(Jac77g)(i^) G (0,oo)}, when Kc < 1, becomes: 

Sa = {^.eM2:(7 + rt,,(^))' + (r;,,(^))'<(i-«.)2} 
= {ueR': {rig{i.),rig{,.)) e ^((-7, 0), i-k,)}. 

We are interested in determining the global expected number of minimum images in D, namely, we shall compute 
(|17p for microlensing. By (J17p . the global expected number of minimum images would then be obtained by averaging 
over the possible uniform star positions (i.e., the randomness arising from the lensing map), followed by averaging 
the source position over a compact set with positive measure S, and then "averaging" over the possible sets in {©}, 
where each has the same area. 

We obtain: 



Theorem 8. Let _D C M^ be compact and suppose that Kc < I. Then the global mean number of minimum images 
lying inside D is given by: 



13 



K, fJ.D,ea f (1 - Kc)^ - {j + zf -w^ [^ ^ Hi{z,w) ^ H3{z,w) 



i.[iV,_(Ay;6)],e> - ^"^J^^-jiri^i^ 



1+ ^^ ' ^ + 



9^ 



dzdw + 0{g ^), 

(18) 
in the large g limit, where 

B = B((-7,0),l - Kc), MD,eo = jg-T: //3(2:,w) = ii'2(2,w;ao), ao se — / |x|^da;, 

and iJi anrf i?2 are defined in equations (Q) and Ol). Recall that \&\ = \&q\ by definition of the global expectation. 

Note that E[Ng„iin{D, Y\ &)\sq\ depends on the common mass m smce H2 is a function of 771. Also the integral 
on the R.H.S. of equation(|18p remains bounded independently of D although oq can be arbitrarily large. 

Proof: By definition, 

E[Ng^^,n{D,Y-&)]{e^ = j^ I E[det{JsiCijg)(x)lg^{x)] dx. 

Ifc'ol Jd 

Using Theorem [TJ the integrand in the R.H.S. of the above equation, which is bounded, can be written in the 
microlensing case as: 

E [det (Jac r,,) {x) Ig^] = E [((1 - K,f - (7 + n,,(x))2 - iTlj\x)) Ig^ (x)] 

{{1- Kef - [j + zf -w^) fr' (^),r- (x){z,w)dzdw 



B 



27r Jb [kI + 



- (7 + z)^-w^ 



9 9^ 



dzdw + 0[g ), 



■z2+w2):V2 
for large g. Integrating over D then yields the desired result. D 

A physical assumption often used by astronomers is the following: On the macro- lensing scale (e.g., galaxy scale), 
the macro-lens potential is averaged over the stars' positions and the lensing map is modeled by the transformation 
•qraacro : D ^ R^ givcn by 

rimacro{xi,X2) = ((1 - Ktot + j)^!^ (1 " l^tot - l)x2), 

where Ktot = k, + k^. In this case, choosing 60 — VmacroiD) gives 

_ 1^1 _ 1 

^^■^° " |6or 1(1 - «:*o02 + 7^1 - ^'""^'^'" 

which we call the macro-average total magnification over D. For the case when 60 = ilmacroiD) , we define: 



Recall that the stars positions are uniformly distributed in the disc B{0, R) of radius R = ^Jragj k,^. 



Corollary 9. Let D^ be a closed disc of radius r . Assume that 60 = VmacroiDr). The global mean number of 
minimum images is: 



(19) 



E[Nmin] = limlimE[Ng^„iin{Dr)] = :^^j^^ 7-^2 721 / ^J; , _^y , ^^^2\R/?. dzdw, 



and the global mean number of images and the global mean number of saddle images grow as order 0{g) for large g, 
with self-evident definitions. 
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Proof: We have fiD,eo = fJ-macro- The functions Hi and H2 defined by equation ([3]) are bounded over R^ and R^ x [0, r] 
respectively for a fixed r. Thus, the function H^ is bounded, independently of g. Therefore, the integrand of the 
integral term in equation (|18p is bounded when we substitute Dr for D. We can then take the limit as g — > 00 in 
equation (jlSp with D replaced by Dr and apply the Dominated Convergence Theorem to obtain 



lim E[Ng^,mn{Dr)] = 



2^1(1 - Ktotr ~ 7^1 



a-^cr-h+z)' 



w 



(«? 



2)3/2 



dzdw. 



which is independent of Dr 
follows from Corollary [6l 



Taking the limit as r ^ (X) yields the first result of the corollary. The second one then 

D 



Remark: 

1. For physically relevant values of k*, Kc and 7, the R.H.S. of (fT9|) is finite. However, the global mean is divergent 

for (1-Ktot)^ =7^- 

2. This corollary agrees with the findings by Wambsganss, Witt, and Schneider 'W\ and Granot, Schrecter, and 
Wambsganss jlOi]. The integral in equation p9|) was resolved in closed form using elliptic integrals (see Ref. 
0)- 

In the table below, we illustrate E[Nmin\ for Ktot = Kc + k* = 0.450 and varying 7 and k,. For all the values of k, 
and 7 shown, we have the "macro-minimum" situation of 1 — Ktot + 7 > and 1 — Ktot — 7 > 0. Note that the number 
of images increases significantly when 7 gets close to 1 — Ktot- In the column for k* — 0.45 = Ktot, note how E[Njnin] 
increases to 6.108 as 7 increases. 





E[Nrmn]-- 


Global Mean Number of Minima 


7\'«* 


0.1 


0.2 


0.3 


0.4 


0.45 


0.1 


1.037 


1.114 


1.219 


1.345 


1.415 


0.2 


1.068 


1.171 


1.230 


1.448 


1.529 


0.3 


1.142 


1.306 


1.488 


1.688 


1.794 


0.4 


1.349 


1.667 


1.985 


2.313 


2.484 


0.5 


2.520 


3.620 


4.621 


5.608 


6.108 



Table I: Global mean number of images that are minima (cf. Corollary [9]) . 

We see from the table that the global mean number of minima is relatively small. For Ktot = 0.450 and admissible 
values of 7, Kc and k», we observe that the global mean number of minima is between 1 and 3. 

Next we present values oi E[Ng^rnin{D)], the global mean number of minima lying inside the disk D = -6(0, 1111.11). 
Here, we &x g = 1000, to = 1 and Ktot ~ 0.450 and we vary the values of k^ and 7. 



E[N 


J,m^n{D)]■. 


Global Mean Numt 


cr of Minima in D 


7\k* 


0.1 


0.2 


0.3 


0.4 


0.45 


0.1 


1.036 


1.113 


1.217 


1.343 


1.412 


0.2 


1.067 


1.170 


1.297 


1.446 


1.526 


0.3 


1.142 


1.305 


1.486 


1.685 


1.791 


0.4 


1.349 


1.666 


1.982 


2.309 


2.479 


0.5 


2.519 


3.617 


4.615 


5.599 


6.096 



Table II: Asymptotics of the global mean number of minima over the disk B{0, 1111.11) (cf. Thcorem[8]). 



The values in Table II differ from the corresponding values in Table I by a relatively small amount. This is expected 
as the difference between E[Ng^rnin{D)] and E[Nmin] has order 0{g~^) (by Theorem [5] and Corollary [5]) . 
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VI. CONCLUSION 

We carried on our development of a mathematical theory for stochastic microlensing, devoting our efforts to the 
study of the expected number of different types of lensed images. This involved an analysis of the random microlensing 
shear. In particular, we obtained the first three leading terms in their asymptotic joint p.d.f.s in the limit of an infinite 
number of stars and showed that each of their marginals converges to a shifted Cauchy distribution, which has a heavy 
tail. The third leading order term of this joint p.d.f. depends on the magnitude of the shear due to the stars, the 
optical depth, and the mean number of stars within a disk centered at the origin. We also obtained the asymptotic 
p.d.f. of the shear's magnitude in the limit of an infinite number of stars. These results on the random shear are 
given for a general point in the lens plane. Using the Kac-Rice formula and Morse theory, we determined the expected 
total number of images and expected number of saddle images for a general random lensing map. By formulating a 
generalization of the Kac-Rice type formula to a family of compact subsets with positive measure of the light source 
plane, we introduced the global expected number of positive parity images for a general random lensing situation. 
Our general results are applicable to a wide range of lensing models. Specifically, we computed the global expected 
number of minimum images for the case of microlensing. This global expectation was also illustrated with explicit 
numerical values. 
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APPENDIX A: PROOF THAT 3„ = 0. 

We consider the integral 

Jr-\x\ ''" 

for any integer n > 2 and < \x\ < R, where 

|a;P + r"^ - R^ 



f{r) — arccos 
Define the variable ■& by 



2\x\r 



Rcosd — rcos/(r) — |a;|, 
with d restricted to the interval [0,7r]. This implies that rsin/(r) = i?sin^. Therefore, 

rexp[i/(r)] = |a;| + Re'^ and rexp[-i/(r)] = |a;| + Re^'^ . 
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We now use these equalities and the change of variable from r to t? to obtain the following: 



^"+1 - 2l 



2i 



R+\x\ 



R-\x\ 

R+\x\ 



R-\x\ 



^i{n+l)f(r) 



-i(n+l)/(r) 



dr 



1 



R\x\ 

R\x\ 
2i 



1 



n+1 



«/('■) 



n+1 



r dr 



{\x\ + i?e-»^)"+' 



(|a;| + Re^^f^^ 



id d{} 



1-& 



yr [X 



Re-^-^y 



-dd 



R\i 



(2i)2i?"+i J_^ 



„i-d 



^-i§ 



(«+l)^ (l 



R ^ 



n+1 



dd 



AR" 



Lt(n+2)i) 



^,ni} ^ ( ^ 



R 



-(n+1) 



kl y- {-fVik + n + l) \x\^ I"' ^^,(„+fc)^ _ ^,(„+2+fc)«^^^ 



4i?" ^^ /cirfn + 1) R'' J ^ 

k=Q ^ ' ^ 

0. 



Above, we used the inequality < |a;| < i? to obtain the Taylor series expansion. Also, the convergence of the series 
allowed us to reverse the order of the summation and the integral. Finally, the last line follows by direct integration, 
where we used the assumption n + 1 > 2, that is, n > 1. 
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